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Abstract 

. We consider the lattice realization of the Standard Model with 

^> \ an additional Z§ symmetry. Numerical simulations were performed 

on the asymmetric lattice, which corresponds to the finite temper- 
^ . ature theory. Our choice of parameters corresponds to large Higgs 

masses (M# > 90 Gev). The phase diagram was investigated and has 
been found to be different from that of the usual lattice realization 
■ of the Standard Model. It has been found, that the confinement- 

deconfinement phase transition lines for the 577(2) and SU(3) fields 
coincide. The transition line between Higgs and symmetric deconfine- 
<^ • ment parts of the phase diagram and the confinement-deconnnement 

transition line meet in a triple point. The transition between Higgs 
and symmetric parts of the phase diagram corresponds to the finite 
temperature electroweak transition/crossover. We see for the first 
time evidence that Nambu monopoles are condensed at T > T c while 
at T < T c their condensate vanishes. 



Owing to the present triviality bound the perturbation expansion in the 
electroweak sector of the Standard Model is considered to work perfectly 
at the energies of the electroweak scale Mz for Higgs masses up to 1 Tev 
[T]. However, it was shown [2], that the finite temperature perturbation 
expansion breaks down at the temperatures above the electroweak transi- 
tion/crossover already for Higgs masses above about 60 GeV. Therefore the 
present lower bound on the Higgs mass, M# > 114 GeV, requires the use of 
nonperturbative techniques while investigating electroweak physics at high 
temperature. Thus our understanding of electroweak physics, which is based 
mainly on perturbation theory could be changed drastically at temperatures 
close to and above the electroweak transition/crossover. 

It was shown recently jHl IH IH] that there is a hidden Zq symmetry in 
the Higgs and fermion sectors of the Standard Model. The Standard Model 
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on the lattice can be denned in such a way, that the whole model is Z 6 
invariant. The resulting model has the same perturbation expansion as the 
usual lattice realization of the Standard Model, which does not respect the 
Zq symmetry. On the other hand, it was argued, that nonperturbatively 
those two lattice models may represent different physics due to their different 
symmetry properties. In particular, it was supposed, that these two models 
may describe the physics at temperatures close to the electroweak crossover 
in different ways. 

In this paper we report results of our investigation of the Z§ symmetric 
lattice version of the Standard Model at finite temperature. We considered 
the model in the London limit, i.e. with infinite bare Higgs mass. This 
does not mean, however, that the renormalized mass of the Higgs boson is 
infinite [6J. 

The phase diagram of our lattice model (Fig. f ) differs drastically from 
that of the usual lattice realization of the Standard Model. Namely, in the 
latter only one phase is present, the phase transition lines degenerate and 
become crossover lines J7]. Our lattice model clearly contains three parts. 
The first one (I) is the confinement phase, where the SU(3) fields confine 
quarks, and confinement-like forces are observed between the leptons. In the 
next one (II) there are no confinement-like forces at all but the line-like ob- 
jects which arise in the unitary gauge are found to be condensed. We identify 
these objects with quantum generalization of the well-known classical Nambu 
monopole configurations 8 . The last part of the phase diagram (III) corre- 
sponds to the low temperature physics, where the Higgs field is condensed. 
In this part of the phase diagram Nambu monopoles are not condensed and 
their density is dropping rapidly when moving away from the transition line. 
Further we refer to both parts II and III of the diagram as to phases, taking 
in mind, however, that the transition line between them could be actually a 
crossover line. 

The increase of temperature corresponds to a shift from phase III to 
phase II. The physical temperature is expressed here as T = — where 
Nt is the time extent of the lattice while a is the lattice spacing, which 
depends on the values of the coupling constants. We denote the value of 
temperature at the transition point as T c . Here the monopole condensate 
plays the role of an order parameter. We do not observe a vanishing of the 
Nambu monopole condensate in phase III with increasing lattice size. This 
result leads us to suggest the hypothesis, that those monopoles survive and 
are condensed in the continuum theory at T > T c . This is in accordance 
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with the supposition which was made in the framework of the SU(2) Higgs 
model in |Tj] l . It was argued in [7] that in the usual definition of the 
lattice Standard Model the electroweak transition is actually a crossover at 
the allowed values of the Higgs mass. This has led to the conclusion that 
the baryon asymmetry could not be produced during the electroweak phase 
transition, as was suggested in ^2]- Our investigation shows, however, that 
the vacuum structure below and above the transition is different. Therefore, 
we do not exclude that there is a phase transition of a high order at T = T c , 
where the condensate of electromagnetic monopoles vanishes. Although the 
high order phase transitions are not well understood, they are known to exist 
in some lattice and statistical systems • The situation may also be similar 
to that of the 3D Compact U{1) Lattice Higgs Model [13], where the phase 
boundary consists of a line of first-order phase transitions at small Higgs self- 
coupling, ending at a critical point. The phase boundary then continues as 
a Kertesz line across which thermodynamic quantities are nonsingular. It is 
worth mentioning, that within the 3D SU(2) Higgs model ^5] it was found 
that the Z-vortices percolate at T > T c while at T < T c they do not. For this 
reason in Ref. ^5] this transition was called the "percolation transition" . It 
is also in accordance with our observations, as in 4D Z-vortices are known 
to terminate at Nambu monopoles. 

We consider here the lattice model described in pj. We use asymmetric 
lattices with time extents 2 and 4 and of space sizes from 8 3 up to 24 3 . With 
the definitions of |5] the pure gauge part of the action has the form 

S g = (3 {2(l-iTr[/ p cos0 p ) + 

plaquettes 

+ (1 -cos2# p ) 

+ 6[1 - iReTrr p Tr[/ p exp(i0 p /3)] 

+ 3[1 - §ReTri>xp(-2zfl p /3)] 

+ 3[l-|ReTrr p exp(4^ p /3)]}, (1) 

1 In JI] it was shown that at a certain limit of the coupling constants the SU(2) 
Higgs model becomes identical to the Georgi-Glashow model. Then, 't Hooft-Polyakov 
monopoles were identified with Nambu monopoles. Therefore, condensation of 't Hooft- 
Polyakov monopoles in the symmetric phase of the Georgi-Glashow model means that at 
least in the limit of coupling constants considered the Nambu monopoles are condensed 
in the symmetric phase of the SU (2) Higgs model. 
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where the sum runs over the elementary plaquettes of the lattice, and 

r G SU(3), U G SU(2), e ie G 17(1). (2) 

Each term of the action, Eq. (JT|). corresponds to a parallel transporter along 
the boundary dp of a plaquette p. 

The action for the scalar field is considered in its simplest form [3] in the 
London limit, i.e., in the limit of infinite bare Higgs mass. After fixing the 
unitary gauge we obtain: 

S H = lY\l-Re{Uiy^)\. (3) 

xy 

Here 7 = v 2 , where v is the bare scalar field vacuum average. We consider 
our model in quenched approximation, i.e., we neglect the effect of virtual 
fermion loops. Thus the whole action of the model is S = S g + Sh- 

It is worth mentioning, that the action of the form Q actually appears 
as a low energy approximation of the SU(5) GUT|3]. Therefore, the bare 
coupling (3 (which is the same for all terms of the action) could be considered 
at the GUT scale. Then, owing to the renormalization group equations, the 
renormalized gauge couplings at the electroweak scale M z come close to the 
experimental ones and coincide with them up to a few percents. Actually, 
their change is of logarithmic form and is very slow. The physical scale (i.e. 
the value of the lattice spacing in physical units) should be determined using 
a measurement of the zero-temperature Z-boson mass in lattice units. We did 
not perform extensive calculation of Mz in our model within the considered 
ranges of bare couplings. However, our measurements of correlators in the 
vector channel allow us to evaluate the lattice spacing in the considered region 
of phase III (J3 G (0.6, 0.8), 7 G (1.0, 1.5)) to be (130 ± 40 GeV)" 1 . 

The renormalized a s (which stands for the strong interaction of quarks) 
in our model can be calculated using certain correlators of colored fields and 
should be expressed as a function of bare couplings. In principle, if we start 
from the theory at the GUT scale, such a calculation must give a reasonable 
result. However, due to the technical problems in lattice simulations, we 
expect direct calculations at the Electroweak scale would give unphysically 
small values of a s (Mz)- This means that at the energies of the order of 
100 GeV color fields in our lattice model appear to be suppressed. Thus 
we consider their influence on the Electroweak dynamics only qualitatively. 
However, if such an influence (which is due to the specific Z e invariant terms 
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in the action) is found, it is reasonable to expect that it also should take 
place for the realistic case of unsuppressed colored fields. 

In general all the renormalized gauge couplings should be calculated us- 
ing static quark and lepton potentials. Then the lines of constant physics 
(LCP) in the space of bare parameters (including time extent of the lattice) 
are defined as the lines, where the zero-temperature renormalized couplings 
and the temperature in physical units are constant [7j. In the present pa- 
per we do not consider LCP, which is needed for a proper investigation of 
approaching continuum physics 2 . The consideration of LCP is also necessary 
for determining the correspondence between the phase diagram in the (3 - 7 
plane and the conventional phase diagram in the T — M# plane. However, 
we can understand the emergence of temperature in the phase diagram rep- 
resented in the Fig. 1 using naive expressions for the lattice spacing and the 
Z-boson mass. Namely, in phase III at tree level Mz = g z v/2 in lattice units. 



Mf hys is about 90 GeV. Therefore, say, at (3 = 0.6, 7 = 1.5 the naive tree 
level estimate for the lattice spacing is a ~ (130 GeV) -1 , while at /3 = 0.8, 
7 = 1.0 it is a ~ (180 GeV) -1 . Finally, the temperature is estimated as 



using the lattice renormalization group equations, it shows, that, in general, 
temperature is increased with an increase of /3 and a decrease of 7 and Nj<. 

Physically interesting values of the coupling constants could be evalu- 
ated following the naive estimates considered above. In our model the bare 
electromagnetic charge is e 2 = g 2 sin 2 6 W = the experimental value is 

e ^ z - > ~ y^8- Thus (3 ~ I ~ 1.9. Our estimate for the critical 7 at (3 — 1.9 is 
7 C ~ 0.9. Therefore, the critical temperature could be estimated as T c ~ 150 
GeV. Of course, this is a very rough estimate and it should be improved using 
direct lattice methods 3 . For technical reasons we did not perform extensive 
numerical simulations in the vicinity of /3 ~ 1.9. Instead we investigated the 

2 We must notice here that for the proper investigation of approaching the continuum 
physics along the LCP it could become necessary to extend the space of bare couplings 
and include there different gauge couplings for SU(2), SU(3), U(l) fields as well as the 
bare Higgs boson mass. 

3 This value should be compared with the one calculated within the SU(2) gauge Higgs 
model in 10 j. There for M H = 120 GeV T c was found to be of the order of 210 GeV while 
for M H = 180 GeV T c ~ 250 GeV 





Although the last expression is to be modified 
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region (3 e (0.6, 0.8) of phase III. This means that our results at the present 
moment should be considered as qualitative only. 

The following variables are considered as creating a Z boson and a W 
boson, respectively: 

Z xy = Z% =sin[Arg^ + ^], 
W xy = W» = Uiy™. (4) 

Here, \x represents the direction (xy). 

After fixing the unitary gauge the electromagnetic U(l) symmetry re- 
mains: 

U X y ¥ 9xU X yQyi 

9 X y -> dxy - a y /2 + a x /2, (5) 

where g x = di&g(e ioix/2 , e~ iax/2 ) . 

In the unitary gauge there is also a U(l) lattice gauge field, which is 
defined as 



rll 

The fields A, Z, and W transform as follows: 



A xy = A* = [-AvgU% + e xy ] mod 2vr, (6) 



A xy ^ A xy a y ~\~ ol x ^ 

7 \ 7 
t-'xy ^xyi 

W xy - W xy e~ ia *. (7) 

It should be mentioned that the field A cannot be treated as a usual 
electromagnetic field as the set of variables A, Z, and W do not diagonalize 
the kinetic part of the pure gauge action ((TJ in its naive continuum limit. In 
our lattice model the electromagnetic field A em should be defined as 

A em = A + Z' - 2 sin 2 6 W Z' , (8) 

where Z' = [AigU xy + 6 l , r j / ]mod27r. The naive value of the Weinberg angle 
corresponds to sin 2 9\y = §■ However, the renormalized Weinberg angle is to 
be calculated through the ratio of the lattice masses: cos 9\y = 

In order to evaluate the zero temperature masses of the Z-boson and 
Higgs boson we use the correlators: 

V F-2/r 

((H x H y )} = (H x H y )-(Hf~—^e- M ^, (9) 
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where H x = J2 y \W xy \ 2 or H x = J2 y Z 2 xy . 

The position of the transition lines on the phase diagram of the finite tem- 
perature model almost coincide with that of the lines on the phase diagram 
of the zero temperature model|)]- The only difference is that the transition 
line between phases II and III is shifted to higher values of 7. Our statistics 
does not allow us to perform a precise calculation of the values of M# and 
Mz- Therefore we have made a very rough estimate of their ratio. Namely, 
in phase III of the zero temperature model (on the lattice 16 ) for the con- 
sidered values of the couplings (/3 G (0.5, 0.8), and 7 up to 1.5) our estimate 
is Mh/Mz ~ 2.5 ± 0.5. This is in qualitative agreement with the predictions 
made within the SU(2) Higgs model considered in the London limit jHj. Thus 
in our model the estimate for the Higgs mass could be M H ~ 230 ± 50 GeV. 
However, as it was mentioned above, the lattice spacing in the considered 
part of the phase diagram is estimated to be (130 ± 40GeV) _1 . In lattice 
theory the inverse lattice spacing plays the role of an ultraviolet cutoff. In 
general, quantum field theory does not work at energies higher than the ul- 
traviolet cutoff. Therefore, we feel it necessary to weaken our estimate for 
M H . Namely, we evaluate it as M H > 90 GeV. 

To understand the dynamics of external charged particles, we consider 
the Polyakov lines defined on the asymmetric lattice in the fermion represen- 
tations listed in the table in [5]: 

V^ pt = (ReTrn {xy)el U xy e- ie ^), 
Vf ept = (ReU (xy)el e~ 2 ^y), 
Quarks = (ReU {xy)el T xy U xy e^y), 

^down quarks — (R- e H(a;j/)gZ Txy e 3 V )i 

^upquarks = (ReU {xy)&l T xy e^) . (10) 



Here / denotes a line on the lattice in the time direction, which is closed due 
to the periodic boundary conditions. 

It is found that to the left of the vertical line of the phase diagram all 
Poliakov lines vanish while to the right of this line all of them increase rapidly 
(Fig. 2). 

In order to extract physical information from the SU(3) fields themselves 
in a particulary simple way we use the so-called indirect Maximal Center 
Projection (see, for example, [T5| IT7j). 
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We investigated several types of monopoles. The monopoles, which carry 
information about colored fields are extracted from the composite fields C l 
(for their definition see [5]) constructed of the SU{2) and U{1) fields and of 
the center vortices appearing in the Maximal Center Projection of the color 
group: 

j C i = ^-*d([dC l ]mod2n). (11) 

(Here we used the notations of differential forms on the lattice. For a def- 
inition of those notations see, for example, |18j.) Pure U(l) monopoles, 
corresponding to the second term in are extracted in the same way from 
2$: J20 — 2^* d([d2d]mod2n) . We refer to them as hypercharge monopoles. 

The electromagnetic monopoles must be related to the field A em . How- 
ever, A ern itself is not a usual lattice U(l) field that should be periodic with 
the period 2tt. Instead A em is constructed of the two U{1) variables: A and 
Z' . Therefore, the electromagnetic monopoles should be constructed of ei- 
ther A or Z' fields, or, possibly both of them. Therefore, we denote these 
monopoles by ja = ^*d([dA]mod27i) and jz = ^* d([dZ']mod27i) . 

The worldsheet of the quantum Z-string may be defined as oz = ^{[dZ'} mod27r— 
dZ'}. This is actually a Nielsen-Olesen string embedded in the Standard 
Model. 

On the classical level the singularity of the hypercharge field 29 = [A + 
Z']mod27r is suppressed by the U{\) pure gauge field action. Therefore, one 
would expect that jz = —Ja- This situation corresponds to the appearance 
of the quantum Nambu monopole with the worldline ja- From (JHJ) it follows 
that its magnetic charge is proportional to Aitsiv^Ow as it should [8 . 

In lattice quantum theory, however, the singularities of the hypercharge 
field may appear in the form of the corresponding monopoles. This situation 
corresponds to the appearance of jie- Then in the absence of a Z-string the 
magnetic charge of such configurations appear to be proportional to 2tt. This 
case of the magnetic monopole for even ]a seems to be corresponding to a 
Cho-Maison monopole or dyon [5] 4 . 

Thus, we arrive at the following two possibilities: 

1. In the absence of j'20 the A-monopoles ja represent Nambu monopoles 
with the magnetic charge AitsivPOy/- 

4 Another way to understand the appearance of a Cho-Maison monopole with the mag- 
netic charge 4-7T is to consider the monopole current extracted from the field 9. Then, due 
to the identity A cm = 29 — 2sin 2 ByjZ' the magnetic charge of such a monopole current is 
47r in the absence of the Z-string. 
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2. If j 2 e 7^ the A-monopoles may represent another type of monopoles 
with the magnetic charge 2ir. Such a monopole with even j& corresponds to 
Cho-Maison monopole or dyon. 

The density of the monopoles is defined as follows: 



where L is the lattice size. 

It is found that the densities of the color and hypercharge monopoles 
decrease rapidly to the right of the vertical line in the phase diagram and 
drop to zero soon after the phase transition. This, together with the behav- 
ior of the Polyakov lines allows us to identify the vertical line of the phase 
diagram with the confinement-deconfinement phase transition common for 
i7(l), SU(2), and SU(3) fields. The behavior of such quantities as the over- 
all action and the monopole densities possess hysteresis effects, which allow 
us to suppose that this phase transition is of the first order. The position of 
the phase transition is localized using hysteresis of the action and is defined 
as the point where the mentioned Polyakov lines vanish. 

The monopole density p(A) constructed of ja is found to be nonzero 
for all values of the couplings considered within phase II. As the hyper- 
charge monopoles disappear in this phase, we identify here ja with Nambu 
monopoles. When going to phase III the density decreases and vanishes soon 
after the transition. In order to investigate the condensation of the elec- 
tromagnetic monopoles we use the percolation probability 11(A). It is the 
probability that two infinitely distant points are connected by a monopole 
cluster (for more details of the definition see, for example, fT7j). We found 
that this probability is an order parameter, which feels the transition (Fig. 3). 
We define the position of the transition using maximum of the susceptibility 
X = (H 2 ) — (H) 2 . It coincides with the point, where the percolation prob- 
ability vanishes. At the same time there is no abrupt change of the action 
on the transition line. The correlation lengths extracted from the space-like 
correlators Q on the asymmetric lattices do not increase when approach- 
ing the transition line between the phases II and III. Let us remind here 
again that in the conventional lattice Standard Model the similar transition 
is found to be a crossover [7] . However, in [15] it was found that the Z-strings 
are condensed at high temperatures in the 3D SU{2) Higgs model, while at 
low temperatures they are not. For this reason this transition was called in 
[To] "percolation transition" . We did not investigated in detail the transition 




(12) 
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between phases II and III in our model. Therefore, we cannot draw a defi- 
nite conclusion as to the nature of this transition. However, the percolation 
properties of Nambu monopoles show that the vacuum structure at T > T c 
differs from the vacuum structure at T < T c . Therefore, we do not exclude, 
that this is actually a phase transition of a high order. 

Investigation of different sizes of the lattices (from 8 3 x 2 up to 24 3 x 4) 
shows that neither the density nor the percolation of Nambu monopoles de- 
crease with increasing lattice size for the values of the couplings consid- 
ered (up to (5 — 2.0) in phase II of the model. Therefore we suppose that 
those monopoles survive in the continuum theory like the Abelian projected 
monopoles of pure nonabelian gauge models. 

The properties of quantum Cho-Maison monopoles are sufficiently differ- 
ent from those of the Nambu monopoles as the hypercharge monopole den- 
sity within the physical phases II and III decrease rapidly when moving away 
from the vertical phase transition line. For this reason we do not exclude that 
Cho-Maison monopoles may completely disappear in the continuum theory. 

To conclude, we have considered the Z 6 symmetric lattice version of the 
Standard Model on the lattice at finite temperature. We must mention, that 
the considered lattice model is obtained as a result of several simplifications. 
First of all, we neglect dynamical fermions. Next, we froze radial fluctuations 
of the scalar field. Finally, our choice of gauge couplings corresponds to 
the SU(5) unified theory. As a result, strictly speaking, the model may 
describe color fields quantitatively only at the energies close to the GUT 
scale. However, we expect that even at the Electroweak scale this model 
may describe qualitatively the influence of the emergence of Z 6 symmetry in 
the lattice action on the phase diagram. We must remind here, that without 
Zq symmetry there is no influence of color fields on the Electroweak dynamics 
at all (if one neglects fermion loops). 

We have found that at least the lattice model itself differs from the 
usual realization of the Standard Model on the lattice. Namely, there is a 
confinement-deconfinement phase transition line common to SU(2), SU(3), 
and U(l) fields. This line and the line of the transition which corresponds 
to the finite temperature electroweak transition/crossover, meet together in 
a triple point. However, we do not consider properly the limit of vanishing 
lattice spacing. Therefore at the present moment we do not exclude that 
the mentioned new features of the Z 6 invariant lattice Standard Model may 
disappear in the continuum limit. 

Nambu monopoles appeared to be condensed at T > T c in our lattice 
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model. We suppose that at T > T c they survive in the continuum theory. 
Finally, if those monopoles are indeed present at T > T c it would be natural to 
suppose, that they may appear at low temperatures in the form of ordinary 
particles 5 . It is worth mentioning, that although we started from the Zq 
invariant form of lattice Standard Model, our prediction of the behavior of 
quantum Nambu monopoles could be unchanged even without taking care of 
Zq symmetry. This should be, of course, the subject of another research. 
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Figure 1: The phase diagram of the model in the ((3, 7)-plane for the time 
extent N T = 2. 
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Figure 2: The expectation value of Polyakov lines as a function of (3 for a 
fixed value 7 = 0.5. (1) corresponds to V^ pt ; (2) corresponds to V^ pt ; (3) 
corresponds to "Pq Uarks . 
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Figure 3: The density p(A) of Nambu monopoles and the percolation prob- 
ability H(A) as a function of 7 for a fixed value (3 = 0.6. 
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